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State-dependent quantum electrodynamic corrections are evaluated for the hyperfine splitting of nS states 
for arbitrary principal quantum number n. The calculations comprise both the self-energy and the vacuum- 
polarization correction of order a (Z a) 2 Ef and the recoil correction of order (Z a) 2 (m/M) Ef- Higher- 
order corrections are summarized and partly reevaluated as well. Accurate predictions for hydrogen hyperfine 
splitting intervals of nS states with n — 2, ... ,8 are presented. The results obtained are important due to steady 
progress in hydrogen spectroscopy for transitions involving highly excited S states. 
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I. INTRODUCTION 

Investigations of the hyperfine structure in light hydrogen- 
like atomic systems are interesting for two main reasons. 
First, accurate measurements of the hyperfine splitting (hfs), 
combined with high-precision ab initio calculations, can yield 
fundamental tests of bound-state QED theory. Second, the 
accurate knowledge of the hfs also constitutes a necessary in- 
gredient in the determination of fundamental constants from 
hydrogen and deuterium spectroscopy. The hyperfine com- 
ponents of transitions in hydrogen can be accurately resolved 
at the current level of spectroscopic accuracy, and the knowl- 
edge of the hfs of excited states is therefore necessary for the 
interpretation of the experimental data. 

The ground-state hfs in hydrogen is known with an out- 
standing accuracy (a part in 10 12 ) for over of 3 decades, and 
the value of 



AE 1S = 1420 405 751.768(1) Hz 



(1) 



has been obtained in Ref. 1 1 ] as a conservative average of var- 
ious experimental investigations of comparable accuracy, the 
first of which was reported in Ref. 1 2] . Unfortunately, our the- 
oretical understanding of the ground-state hfs is limited by the 
insufficient knowledge of the nuclear charge and magnetiza- 
tion distributions, whose contribution of about —50 kHz (30 
ppm) cannot be accurately calculated at present. 

One of the possibilities to overcome this difficulty 0| is to 
study the normalized difference of the nS and IS* hfs inter- 
vals, 



AE. 



nS 



AE- 



is 



(2) 



In this combination of energy intervals, the hfs energy shifts 
due to the nuclear charge and magnetization distributions are 
largely eliminated. Indeed, the lowest-order nuclear correc- 
tions to AEis and AE n s scale with the nonrelativistic elec- 
tron density at the position of the nucleus \4> n {r = 0)| 2 which 
is strictly proportional to n~ 3 . The nuclear effects thus do 
not contribute to the difference A n to leading order. Theoret- 



ical investigations show that the specific difference A„ pro- 
vides an opportunity to test the QED theory of bound states 
on a level of about two orders of magnitude better than for the 
ground-state hyperfine interval AE\s alone 1 1]. According to 
widely accepted terminology, the corrections that depend on 
n through \<j> n (r — 0)| 2 only are called "state independent". 
Thus, only state dependent correction should be considered in 
theoretical investigations of the difference A„. 

Accurate experimental results for the difference A2 are 
presently available for the hydrogen, deuterium, and the 3 He 
ion. Notably, recent progress has been achieved for hydrogen 
0] and deuterium |5] via optical spectroscopy, by comparing 
the 15 and 25 hyperfine splittings via a phase-coherent opti- 
cal measurements of the 1S(F = 0) <^ 25 (F = 0) vs. the 
IS(F = 1) 25 (F = 1) transition. The best absolute ac- 
curacy for the difference A2 is, however, still obtained for 
the 3 He ion in a combination of two relatively old measure- 
ments I 



A 2 ( 3 He + ) = 1 189.979 (71) kHz . 



(3) 



While the specific difference of the 25 and 15 hfs inter- 
vals has been a subject of experimental and theoretical 0, @] 
investigations for a long time, the difference A„ for n > 2 
has attracted much less attention up to now. The case n > 2 
is, however, becoming of significant interest nowadays, due 
to steady progress in hydrogen spectroscopy for transitions 
involving highly excited 5 states. Two ongoing experiments 
could be mentioned in this connection, which concern the hy- 
drogen 15 — 35 transition I9t ll0ll and are expected to reach a 
sub-kHz level of accuracy. 

In the present work, we perform a calculation of the lead- 
ing state-dependent self-energy and vacuum-polarization cor- 
rections for an arbitrary nS state. For the case n = 2, we 
reproduce the well-known results by Zwanziger |8]. We also 
generalize the derivation of the leading state-dependent recoil 
correction given by Sternheim 0] for n = 2 to general n. 
Next, we summarize and partly reevaluate the state-dependent 
higher-order correction and present numerical results for the 
difference A„ with n = 2, ... ,8 for hydrogen. 



2 



This paper is organized as follows: Basic quantities are 
introduced in Sec. |H] Third-order state-dependent correc- 
tions are analyzed and summarized in Sec. |lll] Among these, 
self-energy corrections are treated in Sec. IIII Al vacuum- 
polarization corrections in Sec. lHIBl and recoil corrections in 
Sec. IIII CI The current status of higher-order state-dependent 
corrections is discussed in Sec. lIVI Finally, the total theoret- 
ical predictions for the normalized difference of the hfs inter- 
vals A„ in hydrogen are presented in SecfVlfor n = 2, . . . , 8. 



II. GENERAL FORMULAS AND NOTATIONS 



of the nuclear moment with the magnetic field of the moving 
electron, which in turn is proportional to the orbital angular 
momentum L. We have 



Hm s — H s + H D + H L 
H 8 = &*.p*ir), 



H D = 

Hl = 4" 



|e| 3 a ■ r fl ■ f — a ■ jl 
8m irr 3 
e| jl ■ L 



4m irr J 



(10a) 
(10b) 

(10c) 

(10d) 



We are using natural units with h = c = eo = 1. The 
electron charge is denoted by e = — |e| and a — e 2 /(47r). 
The magnetic dipole moment of the nucleus is 



Here, r is the unit vector in the direction of r. For the 
Schrodinger wave function <j> n of an nS state, the expectation 
value of the nonrelativistic Hamiltonian is 



= 9 Mat / . 



(4) 



where g denotes the nuclear g factor, jU y = \e\/(2m p ) is the 
nuclear magneton, and m p is the proton mass. The vector 
potential generated by the nuclear dipole moment is 



(H Ms ) = (H S ) = 1± (<?• ft) |^„(0)| 2 , (11) 
3m 

and the splitting between the ground-state levels with F = 
I + h and F = I — | gives us the Fermi energy 



4 7T T 47T T 



(5) 



The interaction of the bound electron with the dipole nuclear 
magnetic field is given by the Fermi-Breit operator, 



E F = ^<?MAr|</), i= l(0)| 2 (2/+l), (12) 

3m 

where \<j> n =i (0)| 2 = (Z a) 3 m 3 /tt in the non-recoil limit. 



Vhfs = -e a ■ A = 3 

47r r° 



(6) 



The expectation value of the Fermi-Breit operator on Dirac 
point-nucleus wave functions is well-known. We write it as 

1 



. . 3 g m k 
E hfs = a{Za) ^ — j^ n3 (2k+1)(k 2_ 1/4) 

x A(Z a) [F(F + 1) - 1(1 + 1) - j{j + 1)] , (7) 
where A(Z a) is a relativistic factor [A(Z a) = l+0(Z a) 2 ], 



A{Za) = n 3 \k\(2k + 1) 



2k(7 + n r ) - N 
iV 4 7(4 7 2 - 1) 



(8) 



Here, TV 



yjn 2 , + 2n r 7 + k 2 , 



n — \k\, 7 



^/ k 2 — (Za) 2 , n is the principal quantum number of the elec- 
tron, k is its Dirac angular quantum number, j = \k\ — 1/2 
is the total momentum of the electron, and m is the electron 
mass. 

For future reference, we also give the magnetic field corre- 
sponding to the vector potential @, 



B = V x A= -flS 3 (r) 
3 



3(/I • r)r — fi 



An r 3 



(9) 



In the nonrelativistic limit, the hyperfine Hamiltonian H^is 
is given by the sum of two terms, the first of which is pro- 
portional to a ■ B and is denoted here as Hg + Hp, whereas 
the second one (labeled Hl) corresponds to the interaction 



HI. THIRD-ORDER CORRECTIONS 

A. Self-energy 

The leading state-dependent self-energy correction to the 
hyperfine splitting can be conveniently expressed as 



5Af = -(Z a) 2 Ep {off (n, 1) H{Z a y 2 ] + a^(n, 1)} . 

(13) 

Here, (5A^ E is the contribution to the normalized difference 
A„ due to self-energy effects, where A„ is defined according 
to Eq. (|2j. In general, we will denote various contributions 
to A„ by the symbol SA n with appropriate superscripts. The 
coefficients af^(n, 1) are understood as originating from the 

difference aff(n,l) = af®(nS) - af/(15), with af^{nS) 
being the corresponding coefficient for the nS state. As usual, 
the first index of a?- counts the power of Za, and the second 
one indicates the power of the logarithm ln[(Za) -2 ]. 

The self-energy correction dl3l l consists of two parts in- 
duced by the low-energy and the high-energy virtual pho- 
tons flUl . The low-energy part can be immediately obtained 
by generalizing formulas given in Refs. fill fl2l Il3ll . The 
corresponding contribution expressed in units of a(Z a) 2 /ir 
reads: 



C L 



2 (Za) 2 
x In 



(Z a) 2 m 



!-I + ^ +7 + *(„)-ln(„) 



+ N(nS)-N(lS). (14) 
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TABLE I: Numerical values of the quantity N(nS). 



N(nS) 



1 17.855 67203(1) 

2 12.032 14158(1) 

3 10.449 809(1) 

4 9.722413(1) 

5 9.304114(1) 

6 9.031832(1) 

7 8.840 123(1) 

8 8.697 639(1) 



Here, N(nS) is a delta correction to the Bethe logarithm, 
whose numerical values are given in Table|J] 

We now turn to the contribution due to high-energy vir- 
tual photons. Up to relative order a(Z a) 2 , we can use the 
modified Dirac Hamiltonian -ff ra d (for a derivation see, e.g., 
Chap. 7 of 1 14]), which reads 



p-eFx{W 2 )A +/3m + Fi(V 2 )y 



2m \ 



(15) 



This Hamiltonian leads to various self-energy corrections. 
The first of these is an F 2 (0) correction to the effective po- 
tential, evaluated on the relativistic wave functions. It is gen- 
erated by the following term in Eq. fl5i . 



6H = -F 2 (0)—f3a 
2m 



B = -0(Hs 



H D ) , (16) 



where the Schwinger result ^(0) = a/ (2n) has been used, 
Hs and Hr> are given in Eqs. dlQbt and ( II Oct . respectively, 
and is the Dirac 7 matrix in the Dirac representation. The 
corresponding relative correction to the Fermi energy H2i is 



a_(^\p{H s +H D )\i>) 
2tt 



>\H S \> 



(17) 



Here, ip is the fully relativistic (Dirac) hydrogen wave func- 
tion expanded in powers of Z a, whereas <f> is the nonrelativis- 
tic (Schrodinger-Pauli) counterpart. Under the replacement 
-0 — > 0, Eq. Mil simply gives the leading term a/(2ir). The 
numerator of Eq. d!7i diverges in relative order (Z a) 2 when 
evaluated on an nS state. A finite result is obtained, however, 
when the weighted (or normalized) difference of matrix ele- 
ments is considered. We define the normalized difference for 
the general operator A as 



((A)) = n 3 (nS\A\nS) - (1S\A\1S) 



(18) 



The correction d!7i leads to the following contribution to the 
normalized difference (|2j of hfs intervals, 



Ci 



5A SE - 

_ n 

E F 



a ((4>\0(Hs + H D )\1>)) 
2tt (<I>\Hs\4>) 



(19) 



The second correction (C 2 ) is an F 2 correction to the effec- 
tive potential dl6> . i.e. 



-Fi(O)—0V 2 ff-B, 
2m 



(20) 



to be evaluated on the nonrelativistic wave functions. For the 
third correction C3, we have to evaluate an F{ correction to 
the effective potential dlObk the relevant Hamiltonian can be 
expressed as F{(0) V 2 Hs- The forth correction is a second- 
order correction due to an effective one-loop Lamb-shift po- 
tential, which can be expressed as 



AV = a (Za) 



4 , / m \ 10 



a 
3^ 



Za) 





5" 








+ 6 



(21) 



Here, e is a noncovariant low-energy photon cut-off and V 
denotes the Coulomb potential V — —Za/r. Finally, the 
fifth correction is a second-order contribution due to negative- 
energy states and is induced by the relativistic hyperfine po- 
tential Vhfs as given in Eq. (|6} and the term 



F 2 (Q)—^-E 
2 m 



(22) 



from the modified Dirac Hamiltonian dl5> . where E is the 
electric field generated by the Coulomb potential. From the 
r-scaling of the two involved Hamiltonians, it is clear that the 
resulting operator has to be proportional to 1/r 4 . The prefac- 
tor can be obtained using Dirac algebra and considering the 
fact that the main contribution comes from negative-energy 
states with an energy sa — m. 

The high-energy corrections discussed so far are explicitly 
given by 



Cl =7T- 



Co = 



a U\(3{H S + H D )\^)) 



C, =— 



Ca = 



2tt 


(0 


\H S \<t>) 


a 




» 


12 7T 


(£v 




a 


V2e 






) + 24 



rV 



2a 


HZ 

0? V 


4 




2 m 3 r 4 / 


» ( 





rV 



~^? v {E — H)' 



',V 



rV 



(23a) 
(23b) 

(23c) 

(23d) 

(23e) 



Here, we reemphasize that \ip) is the relativistic wave func- 
tion, \cj)) is the nonrelativistic wave function, and all matrix 
elements (A), by default, are understood in terms of the non- 
relativistic wave function. 

The results for the normalized S'-state difference, expressed 
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TABLE II: Numerical values of the nonlogarithmic self-energy co- 
efficient for the normalized difference [a|^(n, 1)] and for the single 



„SE 
"20 


(nS)] in the range n 


= 1,...,8. 


n 


"20 ( n , J-J 


"20 


i 
1 




1 1 ALL j jo I J(t) 


2 


-5.221 233 33(1) 


11.901 105 41(1) 


3 


-6.705 291(1) 


10.417 048(1) 


4 


-7.402 951(1) 


9.719 388(1) 


5 


-7.809 635(1) 


9.312 703(1) 


6 


-8.076 773(1) 


9.045 565(1) 


7 


-8.266 081(1) 


8.856 258(1) 


8 


-8.407 461(1) 


8.714 878(1) 



are given in Table H7I for n — 1, . . . , 8. In the table, we also 
list the values of a§^{nS) obtained with the help of an im- 
proved 15 numerical value, which we give here for reference 
purposes, 

aff (15) = 17.122 338 75(1). (28) 

This result can be immediately obtained according to the im- 
proved numerical evaluation of the low-energy part as de- 
scribed in Ref. 1 13], which contains a correction to the Bethe 
logarithm induced by a Dirac-delta local potential (see also 
the entries in the forth column of Table II of Ref. IU7I0 . 



in units of a(Z a) 2 /tt, are: 

Ci 19 5 49 



c 2 

C 3 
%{Za) 2 

c 5 

s (Za) 2 



- 1 



19 


5 


48 


+ 8- 






i 


\n 






I 




8 






"(s 


3 






1 


X 


1 — 




n 



_ _ _ [ 7 + *(n)-In(n)l . 



2 /m 
-In — 

3 V2e 



11 

36 



(24a) 
(24b) 

(24c) 



(24d) 



2 1 1 

3 + 2^ + 6^ +7 + ~ ln(n) ' 



(24e) 

Adding all the contributions together, we obtain the follow- 
ing result for the self-energy correction JT51 . 



z(Za) 2 



(25) 



Of course, the dependence on the noncovariant photon energy 
cutoff e disappears in the final answer. The result for the log- 
arithmic term is 1 16] 



aif(M) = ~ 



-- - + _ +7 +tf n -Inn 

4 n 4n z 



(26) 



For the nonlogarithmic term a|o (n, 1), we obtain the general 
result 



<(n,l)= N(nS)-N(lS) 

— ^ + — [7 + *(n)-ln(n)] 



71 79 



48 72 n 144 n 2 36 
3 1 1 



3 ln(2) 



--- + — + 7 + *(n) - In(n) 
4 n 4n z 



In the particular case n = 2, we reproduce the known value for 
this coefficient 1 15]. Explicit numerical results for a^o^n, 1) 



B. Vacuum polarization 

The leading state-dependent vacuum-polarization correc- 
tion to the hyperfine splitting can be conveniently expressed 
as 



6 A™ = -(ZaYE F a^(n,l). 

7T 



(29) 



The correction <5A^ P consists of two parts 0], with the first 
one given by a matrix element of the radiatively corrected ex- 
ternal magnetic field and the other by a matrix element of 
the vacuum-polarization operator between the wave functions 
corrected by the presence of the external magnetic field. 

We start with the first part. To the leading order, the radia- 
tively corrected magnetic interaction (magnetic loop) is well- 
known to be 

VvP,mag(f) = Vhfe(r) 



(30) 



We recall that the matrix element of Vm s between the Dirac 
wave functions is, for nS states, 



(n | Vhf s | n) 



E F 



m 2 (Z a) 3 J 



dr g n {r) / n (r) , (31) 



where g n and f n are the upper and the lower radial component 
of the Dirac wave function, respectively. We thus immediately 
have that 

^VP.mag = ( n |^ | n ) 



e f 2a r°° yW^i f j_ 

m 2 {Zaf 3vr J 1 1 t 2 \ + 2t 2 

x / dr(l + 2mrt)e- 2mrt g n {r)f n {r). 
Jo 



(32) 



(27) 

To the leading order in Z a for an nS state, 



9n{r)= 2 -( P -\ e-^CUr-^), 03) 



.3/2 



1 (W 



5 



and 



Zm ar 



(34) 



where (3 — Z am, and £} l _ 1 are generalized Laguerre poly- 
nomials. Performing the integration over r in Eq. (I32> with 
help of entry (2.19.14.6) in Vol. 2 of Ref. 0, expanding the 
result in Z a, and integrating over t, we obtain 



8EY<™z = ^~(Za) 
n 7r 



8 15 



The corresponding contribution to A n is 

<5A VP,ma g = ^ (Za)2i?F _| , , 



(35) 

(36) 



The second vacuum-polarization contribution is given by 
the second-order correction, 



SE, 



VP, el 



Vhf, 



1 



(E - H)> 



VP 



(37) 



where Vvp is the vacuum-polarization potential. Due to 
spherical symmetry of Vyp, only the nS intermediate states 
contribute in the above expression. To the leading order, we 
have 



Vy P {r) 



4 a(Z a) 
15 m? 



S(r) 



(38) 



and we can replace Vhta — > Hs, with Hs being given in 
Eq. dlObt . The second-order matrix element d37t diverges for 
nS states. It is, however, finite for the normalized difference, 
with the result 



5A 



VP.cl 



o 7ra 
~ ~15 F m 5 (Za) 2 



1 



(E-H 



(39) 



Using the formulas from Ref. 1 17] for the matrix element, we 
arrive at 



SAZ P - cl =-(Za) 2 E F 



7T 

X I - 



15 



1 h 7 + *(ri) - ln(n) 

n 



(40) 



Finally, the total result for the vacuum-polarization correc- 
tion [Eq. d29l l reads 



«20 P ("-: 1) 



15 



(41) 



in agreement with Ref. Jl6ll . 



C. Recoil corrections 

The leading-order state-dependent recoil correction can be 
parameterized as 

^H^f^.ctM), (42) 



where M is the mass of the nucleus. The general expression 
for this correction was derived by Sternheim jj]. It reads 



SE* EG 
where 



H 



rr(l) _ 
M - 



H 



(2) 
M 



8m 2 

e 



2iJ 



V • E 



(i) 

M 



rr(2) 
7 W M 



—3 p-SA, 

8to j m 



2m 



a-B, 



(43) 

(44a) 
(44b) 



H (3) _ 



p2 B-B. p2 



Am 2 



Am 2 



2m a ' i 

— ExA+ — (SEx p-px SE) 
2m Am y ' 

(px A — Ax p),p 2 I . (44c) 



8mM 



Here, A is given in Eq. 0, (5E 1 is the electric field induced by 
the scalar potential of a moving magnetic dipole SV, 



8V = ~ — 
An 



„ Ze -A p -1 

— I X TT ' V- 



2M 



M 



and 5 A is the vector potential of the moving nucleus, 



1 Za I _ r r 



(45) 



(46) 



The matrix elements in Eq. J43l > diverge for nS states, but they 
yield a finite result for the normalized difference A„, which 
reads 



,REC 
i 20 



1 - 



1 



7,/ 



17 



1-^7 
2r] 



9 


1 


14n 


28n 


11 


1 


~12 ' 


f 2^ 



12n- 



+ 7 + *(n) - ln(n) 



(47) 



where 77 = <?M/ (Zm p ) and m p is the proton mass. For the 
particular case n = 2, our result is in agreement with the one 
originally obtained by Sternheim yj. 



D. Summary of the theory up to third order 

To the leading order in the parameters a, Z a, and m/M, 
the normalized difference of the hyperfine-structure nS inter- 
vals A n is given by the sum of the relativistic (Breit), self- 
energy, vacuum-polarization, and recoil corrections: 

A n = (Za) 2 E F \a^(n,l) + -\4f(n,l) \n[(Z a)~ 2 } 
+ 4*(n, 1) + 4 p (n, 1)] + g 4f c (n, 1)} , (48) 



6 



where the Fermi energy Ep is defined as the splitting between 
the ground-state levels with the atomic angular momentum 
F = I + 1/2 and F = 1-1/2 calculated within the non- 
relativistic approximation and is given by 



4 3 m 2 /! 21+1 

E F = -a{Z af —— 

3 m p fiN 21 



1 + 



M 



(49) 



with the nuclear magnetic moment u = g fx^ I. Notice that 
this expression follows from Eq. (I12> after restoring the cor- 
rect reduced-mass dependence. 

For the particular (and the most important) case n = 2, 
the coefficients in Eq. J48l > were obtained long ago 
The full n dependence of the coefficients aff and a^o was 
reported in Ref. 1 16]. In the present investigation, we have 
derived the results for all coefficients in Eq. J48i for general 
n. The self-energy, vacuum-polarization, and recoil correc- 
tion are given by Eqs. ( I27l i. J41I . and j47l >. respectively. The 
remaining second-order Breit contribution to A„ is given by 



a 2oK !) = 



3 

2n 



11 

6 n 2 



(50) 



IV. HIGHER-ORDER CORRECTIONS 



According to Ref. fl^l . analogous considerations are valid 
also for the radiative-recoil correction, and hence 



^rad-n* = £ {Za f 111 Ep ^[{Zo)^} 



16 

t 



m 
M 

3 1 

4 ~ ~ 



n An 2 y ' w 



(53) 



We now turn our attention to the state-dependent recoil cor- 
rection to order (m/M) (Z a) 3 Ep, which we evaluate in the 
logarithmic approximation. We have identified two such con- 
tributions. The first one can be obtained as a second-order 
perturbation correction induced by two effective local poten- 
tials, the first one being Hg [Eq.l llObH and the second one 
corresponding to the logarithmic recoil correction to the Lamb 
shift to order (Z a) 5 m 2 /M. The result is 



SA 



HREC.a 



{Zaf 

7T 

4 



-Epln(Za) 

1 - - +7+ - ln(n) 
n 



(54) 



Higher-order QED and nuclear corrections to the difference 
Aq were extensivel y in vestigated during the last years 
Eot l2ll I22I I23I I24I l25ll . The general n dependence of the 
difference A„ received significantly less attention up to now. 
In this section, we would like to summarize the results for 
higher-order corrections and reevaluate some of them. 

The higher-order relativistic (Breit) corrections are imme- 
diately obtained by expanding the general formula 



E F v ; V36 Sn 



67 



55 



36 n 2 12 n 3 
245 721 1195 



21 

8^ 



216 ' 48 n 432 n 2 
33 147 163 



144 n 3 



16 n 4 16 n 5 



(51) 



where the sixth-order contribution is included for complete- 
ness. 

The state-dependent two-loop correction to order a 2 (Z a) 2 
was found in Ref. 1 16] in the logarithmic approximation. This 
result can be easily derived if we observe that the leading one- 
loop aio correction for the ground-state hfs is generated by 
an effective magnetic form-factor correction [Eq. J23aH to the 
Hamiltonian JlObi . We thus employ JlObl as an input for a 
Dirac-delta correction to the Bethe logarithm and obtain the 
result 



c A two-loop 

3 

4 ~ 



(ZafEp ln[(Zay 2 } 

+ -^r + 7 + *(n) - ln(n) 
4n z 



(52) 



in agreement with Ref. lll6ll . 



This expression generalizes the result for the difference A 2 
reported in Ref. 12311 . The second contribution (absent in 
Ref. [23i ) is obtained as a second-order perturbation induced 
by the operator Hg and by the operator responsible for the 
nonlogarithmic recoil correction to the Lamb shift to order 
(Za) 5 m 2 /M. The logarithm of Z a then arises from the 
second term of the Z a expansion of the electron propagator 
after an integration over the logarithmic region |26]. The re- 
sult reads 



SA 



HREC.b 



(Za) 



M 



2,s 



E F ln(Za) 
.i + i-+ 7 + ¥(n)-ln(n) 



(55) 



We note that this contribution, unlike Eq. d54b . is finite for 
single nS states. For 15 state, the constant in Eq. d55t turns 
into (124/9 + 28/3 In 2), which coincides with a part of the 
complete 15* result obtained by Kinoshita l27il (2 Cs in his 
notation). Our result for the logarithmic part of the fourth- 
order recoil correction is the sum of Eqs. (I54i and d55i . 



SA 



HREC 



(Za) 3 



— E F In(Za) 

3 
4 



+ 1- +7 + *(n) -ln(n) 
in 



(56) 



We do not have a proof that this result is complete. 

Some incomplete results for the fourth-order one-loop self- 
energy and vacuum-polarization corrections were obtained in 
Ref. Ilal . With misprints being corrected in 1 22], these cor- 



7 



rections read, respectively, 



5A^ = a(Za) 6 E F 
9 



11 

20 



1 

n 



20 n 



621 n 2 - 1 
" 320 n 2 

2 + 7 + *(n) 



191 
~L6~ 

ln(n) 



- 5 In 2 

(57a) 



HVP 



a {Zaf E F 



55 
26 



1 

n 



81 
26 n 2 



13 
"24 

+ 7 + - ln(n) 



(57b) 



It should be noted that the one-loop self-energy correction 
yields the largest contribution among all fourth-order correc- 
tions mentioned so far and the incompleteness of the result 
( 15 7 at provides the dominant theoretical uncertainty for A„. 
For the particular case n = 2, this correction was evalu- 
ated numerically to all orders in Z a in Refs. The 
deviation of the contribution (I57a> from the all-order result 
was found to be on the level of 20%. The evaluation of the 
complete result for the fourth-order vacuum-polarization cor- 
rection is a much simpler task than for the self-energy. It 
can be solved either analytically, as was done for n = 2 in 
Ref. |22l HH, or (which is much easier) numerically, as was 
done for n = 2 in Ref. 12511 . However, in view of the absence 
of complete results for the self-energy correction, we do not 
pursue the matter any further in the current investigation. 

The nuclear-structure correction was found in Refs. flfilEcIl 
to be 



SA 



Nucl 



-(ZaYAE^ 



Nucl 



5 1 9 

4 n An 2 



+ - ln(ra) 



+ -{Zaf 



n — 1 / Rm \ 2 n 2 — 1 



Re J 4n 2 



7 + *(n) - ln(n) 



(m Re) Ee , (58) 



where Re and Rm are the electric and the magnetic charge 
radii, respectively, and AE^" ' is the nuclear correction for 
the ground-state hfs. 



V. THEORETICAL RESULTS FOR A„ 

In this section, we collect all theoretical contributions avail- 
able to the normalized difference of nS states A„ [Eq. (0]. 
Numerical results for individual contributions and the total 
theoretical values of A n in hydrogen are listed in Table |in| 
for principal quantum numbers n = 1, . . . , 8. The second- 
and third-order corrections summarized by Eq. d48i are given 
in the first five rows of this Table. Forth-order QED correc- 
tions discussed in Sec. lIVI are tabulated in the next seven rows, 
and the nuclear-structure correction completes the analysis. 
Parameters of the proton used for calculating numerical data 
in Table [HI] agree with those from Table 8 of Ref. id . The 



nuclear-structure correction for th e gr ound-state hfs that en- 
ters Eq. d58> was taken from Ref. 12211 . where it was obtained 
by subtracting all known QED corrections from the experi- 
mental result for the ground-state hfs 0. Its numerical value 
is -46 kHz. 

We already mentioned above that in the particular case 
n = 2, there are complete all-order results available for the 
<5A^ SE and <5A^ VP corrections. We thus employ the numer- 
ical values for the self-energy and vacuum-polarization re- 
mainder functions for the difference A2 as given in Ref. 1 25], 
as well as the uncertainty estimates given in the cited refer- 
ence. The corresponding entries in the table are marked with 
the asterisk. For n > 2, we use the formulas d57al and (I57b> 
and ascribe the 50% uncertainty to them. The error estimates 
for the other forth-order corrections are as follows: for the 
two-loop and the radiative recoil corrections, we assume the 
uncertainty to be a half the numerical value of the logarith- 
mic terms, while for the recoil correction we use 100% of the 
correction given by Eq. d56i . 

The two last rows of TablelrTTlare reserved for the total the- 
oretical predictions for the normalized difference A„ and for 
the complete values of the hfs frequency of excited hydrogenic 
nS states. The latter are obtained by combining the highly ac- 
curate experimental value of the ground-state hfs interval Q 
and the theoretical prediction for A„ given in the previous row 
of the table. 

For the case n = 2, our evaluation differs from the previ- 
ous investigation of the difference A2 presented in Ref. ll22ll 
in two ways: (i) we employ the latest numerical results for 
the self-energy remainder from Ref. Il25ll and the error esti- 
mate from this reference and (ii) we also have found an ad- 
ditional (numerically small) higher-order logarithmic recoil 
contribution ( I55> . Despite the small change of the theoret- 
ical prediction, our final result for the hfs frequency of the 
25 state still deviates by 1.4 a from the experimental result 
E 2S = 177 566 860(16) Hz H. We mention also a similar 
(1.8 cr) deviation of the theoretical value of A2 for the 3 He 
ion from the experimental result Q observed in Ref. 1 25 ] . 



VI. CONCLUSION 

The normalized difference of the hfs intervals A2 = 
8 AE2S — AEis has been a subject for both theoretical and 
experimental investigations since a long time. In this paper, 
we have presented calculations that generalize the previous 
studies of A„ = n 3 AE n s — AE\s to general n. Our re- 
sults are complete through third order in the parameters a, 
Z a, and m/M; an estimation of the fourth-order corrections 
is also supplied. 

The dominant source of the present theoretical uncertainty 
for the difference A n comes from the higher-order one-loop 
self-energy correction. Further improvement of the theory can 
be achieved by a numerical all-order (in Z a) evaluation of 
this correction. Such a calculation has been carried out for the 
difference A2 in Refs. Il2j|.|25j based on a method developed 
by a number of authors I2S& l30l 13111 and seems feasible for 
higher values of n as well. It should be noted that the results 



8 



for hydrogen reported in Refs. I2ll l251 involved an extrapola- 
tion of numerical data obtained for Z > 5 towards Z = 1. It 
would clearly be preferable to perform a direct numerical cal- 
culation of the higher-order self-energy correction for Z = 1, 
as it was done for the Lamb shift in Refs. 13 2l 13311 . This project 
is underway. 
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TABLE III: Individual contributions to the normalized difference A n of hfs frequencies, and absolute values of the hyperfine splitting 
frequencies of excited 5* states in hydrogen. For the entries marked with an asterisk (*), we employ the numerical results for the self-energy 
and vacuum-polarization remainder functions as reported in Ref. 1 25] instead of the analytic expressions given in Eqs. I57ai and <57bl used in 
other cases. The absolute values for the hfs frequencies of excited states are obtained with the help of 15 experimental result in Eq. Q as a 
reference. Units are Hz. 



Effect 


2S 


3S 


iS 


5S 


6S 


7S 


8S 


(Za) 2 


47 222.0 


47 571.8 


44 860.9 


42 310.9 


40 226.1 


38 548.6 


37 187.3 


a (Z a) 2 (SE) 


1 936.0 


2718.6 


3 134.2 


3 390.9 


3 564.9 


3 690.4 


3 785.3 


a {Zaf (VP) 


-58.0 


-79.2 


-90.1 


-96.8 


-101.3 


-104.5 


-106.9 


(Za) 2 (m/M) 


-162.9 


-210.3 


-232.6 


-245.6 


-254.0 


-260.0 


-264.4 


Sum of 3 rd order 


48 937.1 


50 000.9 


47 672.4 


45 359.4 


43 435.7 


41 874.5 


40 601.3 


(Za) 4 


5.6 


5.6 


5.2 


4.9 


4.6 


4.4 


4.2 


a 2 (Za) 2 


3.3(1.7) 


4.5(2.3) 


5.1(2.6) 


5.5(2.8) 


5.8(2.9) 


6.0(3.0) 


6.1(3.1) 


a (Za) 2 (m/M) 


-3.1(1.6) 


-4.2(2.1) 


-4.8(2.4) 


-5.2(2.6) 


-5.4(2.7) 


-5.6(2.8) 


-5.7(2.9) 


a (Z a) 3 (SE) 


9.7(5)* 


15.8(7.9) 


19.1(9.6) 


21.2(10.6) 


22.7(11.3) 


23.7(11.9) 


24.5(12.3) 


a (Za) 3 (VP) 


3.0* 


3.7(1.9) 


3.8(1.9) 


3.7(1.9) 


3.7(1.9) 


3.7(1.8) 


3.7(1.8) 


(Za) 3 (m/M) 


0.3(3) 


0.4(4) 


0.4(4) 


0.5(5) 


0.5(5) 


0.5(5) 


0.5(5) 


Sum of 4 th order 


18.7(2.3) 


25.8(8.7) 


28.8(10.4) 


30.6(11.4) 


31.8(12.2) 


32.7(12.7) 


33.3(13.1) 


Nucl 


-1.8 


-1.8 


-1.7 


-1.6 


-1.5 


-1.5 


-1.4 


Total A„ 


48 954.0(2.3) 


50 024.9(8.7) 


47 699.5(10.4) 


45 388.4(11.4) 


43 466.0(12.2) 


41 905.7(12.7) 


40 633.2(13.1) 


HFS freq. 


177 556 838.2(3) 


52 609 473.2(3) 


22 194 585.2(2) 


11 363 609.1(1) 


6 576 153.79(6) 


4 141 246.81(4) 


2 774 309.35(3) 
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